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Abstract 

Thin film alloys with perpendicular anisotropy were studied using Lorentz Transmission Elec- 
tron Microscopy (LTEM). This work focuses on the configuration of domain walls and demon- 
strates the suitability and accuracy of LTEM for the magnetic characterisation of perpendicular 
magnetic anisotropy materials. Thin films of chemically ordered (LIq) FePd alloys were inves- 
tigated by micro-magnetic modelling and LTEM phase retrieval approach. The different com- 
ponents of magnetization described by the modeling were studied on experimental images and 
confirmed by LTEM contrast simulation. Furthermore, quantitative measurements of magnetic 
induction inside the domain walls were made by using an original method to separate the elec- 
trical and magnetical contributions to the phase information. Irregularities were also observed 
along the domain walls which could play a major role during the magnetization processes. 
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Introduction 



Lorentz microscopy is used to reveal the magnetic configuration of magnetic domains 
in a specimen. This powerful technique offers a sensitivity of a few tens of nm.T (this 
value means that more or less 10 nm of materials with a saturated magnetic induction 
of 1 T is needed to be observed by LTEM) coupled with the spatial resolution of electron 
microscopy (chap. 2 of III]). The basis of LTEM is the deflection of the electron beam 
by the Lorentz force during transmission through a magnetic specimen. For this reason, 
only magnetic induction perpendicular to the electron beam can be detected. For in-plane 
magnetization, a classical LTEM observation (Fresnel and Foucault contrast) reveals black 
and white contrast representative of domain wall localisation or domains orientation ||2(]. 
In the case of thin foils with perpendicular magnetization {i.e. magnetization parallel to 
the electron path), this contrast can be obtained by tilting the sample to cause magnetic 
induction perpendicular to the beam |3, 4, ^. However, in the perpendicular magnetical 
anisotropy (PMA) materials class, the domain walls (namely Bloch walls) exhibit an in- 
plane magnetization component. The specimen could therefore remain untilted, with the 
magnetization of domains parallel to the electron path and all deflections occurring in the 
domain walls. To obtain a LTEM signal localized in the domain walls, they must be larger 
than the spatial resolution of the Lorentz lens or the lens used as an imaging lens (from 
less than 1 nm for a dedicated Lorentz lens to few nm for a traditional objective minilens). 
Moreover, the induction perpendicular to the beam must be strong enough to significantly 
deflect the electron beam (starting with a magnetic sample with a saturated magnetization 
of [Iq.Ms = 1 T, a thickness of pure perpendicular-to-the-beam magnetization of more 
than 10 nanometers is needed to deflect significantly the electrons, about 5 [^irad in that 
example [2]). In this paper, we present the study of FePd alloys, which, when chemically 
ordered, exhibit a strong perpendicular anisotropy leading to an up and down domain 
configuration |^. Domains are separated by Bloch walls with an in-plane induction 
around 10 nm in width jZi]. This is large enough to be detected by Fresnel contrast and 
analysed by mapping methods, as described in the next section. 

Fresnel and Foucault contrasts are both standard imaging modes but neither of them gives 
direct quantitative information about the magnitude of magnetization. Image series are 
needed to quantify the observed deflections. Differential Phase Contrast is a mapping 
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technique that has been implemented initially in a STEM microscope la] and then extended 
to conventional TEM ||90. This technique is an extension of imaging by Foucault contrast 
and has not been used in this work. Using a wave-mechanical approach, the action 
of the magnetic induction on the electron wave is modelled as phase shift. This is the 
well-known Aharonov-Bohm effect [lOO. Although phase information is lost by image 
formation, focal series allow phase retrieval by the solving of the Transport of Intensity 
Equation (TIE) 111]. The TIE equation, 

(1) 



relates the phase shift to image intensity along the optical axis, with A being the wave- 
length of the electrons, lx,y being the intensity distribution, (p being the phase shift of the 
electron wave and V^-y is the in-plane gradient |12l] . Nowadays this technique has been 



widely developped, by means of numerical calculations |13l,ll4|l or precautions in the use 
of mathematical approximation ||15l. The technique has been applied to a wide 

variety of both magnetic |1^ arid non-magnetic materials |20]. 

The variation of intensity along the z axis is related to the gradient of the phase shift. As 
it has been described by Aharonov and Bohm, the phase shift is the resuolt of two contri- 
butions : an electrical one, due to the mean inner potential of the material and a magnetic 
one, due to the magnetic induction distribution in and around the sample, namely 

(\){x,y) = CE J Vmt{x,y)dz-^ ^ A^{x,y,z)dz, (2) 

where Ce is electron energy (E) dependent (Ce = 2ne/AE ■ (E + Eq)/{2E + Eq) with Eq = niQC^ 
the rest energy of electrons), Vint is the mean inner potential of the sample (the electrical 
phase shift origin). Az is the z-component of the magnetic vector potential describing the 
magnetic induction distribution in a plane (for a given z) perpendicular to the optical axis 
(B = V X A). 

Knowledge of the phase shift gives access to quantitative magnetic information if the 
electrical contribution can be removed. Various techniques have previously been used to 
achieve thislzil]. Most of these methods need the acquisition of two different phase in- 
formations (changing the sample orientation, saturation condition or accelerating voltage 
for example). In this work, a novel, easily implemented method, with only one phase cal- 
culation is proposed for thin foil samples prepared with the classical method (mechanical 
polishing and ion milling). 



Calculation 



Consider the definition of the magnetic vector potential, the z component variation in 
the X, y directions can be expressed by : 



dyA^ = B, + d,Ay 



d,A, = -By - d,A,. (3) 
Integrating over all the z-axis directions yields, 

1 ^yAz(x, y, z)dz = I B:,{x, y, z)dz - [A,,.(x, y, 
I dAz{x,y,z)dz= I -By{x,y,z)dz-\Ay{x,y,z)\ 

kJ oo U OO 



1Z=oo 
Z=— oo 



(4) 



Assuming that magnetic vector potential vanishes at the infinity (or is the same on each 
side, far from the sample), one gets 

^kx h{x, y) = - V • (/) (x, y) + CeV • [V{x,y)t{x, y)] (5) 

where is a unit vector of the optical axis. The main effect of the vectorial product on the 
l.h.s. of Equation|5]is the inversion the x and y components of b{x, y) = J^B{x, y, z) dz, the 
integration along the optical axis of the magnetic induction including magnetization in 
the sample and stray fields. V{x, y) is the integration of the mean inner potential over all 
the thickness of the sample. 

It is clear that for a sample of constant thickness and composition, the electrical contri- 
bution disappears, and for a uniform gradient of thickness and uniform composition, the 
electrical contribution is reduced to a constant. Hereafter, we show how this equation 
can be used to remove the electrical contribution from experimental phase measurements 
obtained using conventionnal electron microscope specimens. 



Experiment 

This paper presents results obtained from a sample with PMA configuration, exhibiting 
two different in-plane contributions. Previous LTEM observations have been carried out 
on similar alloys Js] and the authors reported that it was necessary to tilt the sample in 
order to reveal the domain structure. In our case the in-plane components lead to strong 



Fresnel contrast without tilting and we were able to reveal the domain wall configuration 
using reconstruction by TIE solving. The information given by the phase retrieval is thus 
directly related to domain wall information and not to the domain magnetization, parallel 
to the electron beam. Such dependance enables a quantitative exploration of the domain 
walls rather than the whole sample. 

The sample was grown on a MgO (001) substrate by Molecular Beam Epitaxy (MBE) 
according to the following sequence : a thin Fe layer (2 ML) was deposited in order 
to initiate the epitaxial growth on the substrate, al6 nm chemically disordered FePd2 
layer was first codeposited as a magnetic layer having a vanishing anisotropy. Then, a 
FePd ordered (LIq) layer (37 nm) with a strong PMA was added at 350°C. The film was 
then covered with a 1.5 nm Pt capping layer to prevent oxidation. The structural and 
magnetic properties of this type of alloy can be found elsewhere j^. Sample for electron 
microscopy was prepared using a classical method by mechanical polishing and Ar ion 
milling. The preparation was in a plan- view geometry and the angle used for ion milling 
was 6 degrees. The microscope used for the study was an FEI TITAN equipped with a 
dedicated Lorentz lens operating at 300 kV and fitted with a Tridiem GATAN Imaging 
Filter. All images were zero-loss filtered and data processing has been completed by 
scripts we have developped in Digital Micr ogr aph from Gatan. The TIE operation was 
carefully done using the fourier approach|l3l], enhanced by image symetrization|l4l]. 
Micromagnetic simulations were performed with GL_FFT (copyright CNRS, Institut Neel 
Grenoble |j22D). The micromagnetic configuration was used to simulate Fresnel contrast 
and compared to experimental images. 



Results 



The focus here is on micromagnetic simulations carried out for this magnetic stacking 
(see Fig 1). The simulations are based on the temporal integration of the Landau-Lifshitz- 
Gilbert equation using regular prisms^ to discretize the sample. The input parameters 
used for FePd were chosen as follows [6i] : Exchange constant A = 6.9 10"^^ J-m"^, Uniaxial 
Anisotropy K = 1.03 lO^J.m"'^, Saturated Magnetization /ioMs = 1.294 T. The cells are 
0.781 nm X 0.625 nm and infinite along the y direction. The results are presented here in a 
cross-sectional view, with the growth direction along the z axis. 
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The top magnetic layer with a LIq structure can now be observed. This layer exhibits an 
out-of -plane magnetic anisotropy which can be described as a succession of up and down 
domains separated by Bloch walls (mark 1 on Fig[T]) and with a periodicity of 160 nm. 
These domain walls are surrounded by Neel Caps due to the demagnetizing field. The 
bottom magnetic layer is the soft layer, where the magnetization tends to be in-plane. 
This implies a larger Neel Cap below the domain walls (mark 2). In-plane components 
of the Neel Caps are usually considered to be symmetric : in terms of electron deflection 
the two deflections are cancel one another. In this sample the balancing no longer exists. 
Below the domains (mark 3), the magnetization enhancement shown previously (mark 2) 
gives rise to another in-plane component. This component is parallel to the domain walls 
and generates a magnetic deflection. Considering the in-plane induction sensitivity of 
LTEM, three magnetic components considered to be acting on the electron beam in this 
sample : (1) Bloch walls, (2) bottom Neel Cap, (3) soft layer under the domains. In our 
2D case, the second term deflects the electrons along the y direction. The technique is 
not sensitive to that magnetization component because the deflection happens along the 
Bloch wall. It would only be apparent at the point of a sudden change in the wall direction. 

Magnetic Fresnel contrast is observed on the defocus image (Fig |2lB) for an untilted 
sample. This contrast comes directly from the domain wall area (mark 1 on Fig |2lC) 
where the magnetization lies in-plane. No deflection occurs in the domains themselves. 
The typical Fresnel contrast for a Bloch wall is a couple of bright and dark bands 
revealing a continuous variation of in-plane magnetization (see Fig|2lD). As shown by 
the simulation, contrast coming from the soft layer (mark 3) is hidden by the intensity of 
the Bloch wall contrast. 



We performed a phase retrieval process via the solution to the TIE (FigSJ. The method 
employed here was the numerical method described by De Graef et al. in 1^. Over- 
and under-focused images used for the reconstruction are presented in Fig SK-B. The 
calculated phase (Fig IH-C) exhibits two kinds of variation. The first is a continuous 
gradient along the phase image of about 10 mrad (along the diagonal of Fig HlC, from 
dark to bright contrast). Smaller variations are also observed in the same image, which 
are measured around 0.1 mrad. The gradient in the vertical and horizontal direction gives 
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a vectorial distribution which can be displayed as a scattergram (inset of Fig IHC). The 
scattergram is a 2D representation of the distribution of the magnetic induction vectors, 
i.e. a bidimensionnal histogram of the magnetic vectors distribution (see also appendix). 
As deduced from Eq. |5] there is both an electrical and a magnetic contribution to this 
scattergram (see also FigHlE). Stated a constant variation of the thickness profile due to 
the sample preparation (this has been checked using the inelastic electron ratio ||23tl on 
the sample, results not presented here) the electric contribution is reduced to a constant 
and the magnetic information is therefore described by 

(pMag = (pMeasnred - (pEl (6) 

kxBl = V^c.-cste. (7) 

> 

Vcak. describes the vector field deduced from the phase gradients (and presented as a 
scattergram in the inset of FigHlC). The cste represents the vectorial contribution (derived 
from the phase) which has no magnetic origin. 

In the scattergram, the gradient of the magnetic phase leads to magnetic induction 
vectors (with coordinates that are rotated by 90° with the use of kx) while the gradient of 
the electrical phase is the same for each pixel of the phase image. As a result, each pixel 
of the gradient-phase image can be described as a sum of a constant vector due to the 
electrical contribution and a vector due to the magnetic contribution. A straightforward 
way to remove the electrical contribution is to remove the vectorial constant due to the 
electrical contribution (see the scheme in Fig IHE). This is done here by considering an 
isotropic distribution of magnetic moments in the image, which has been checked on 
previous MEM measurement | i24( ]. Up and down domains are interlaced so the domain 
walls orientations are pointing in all in-plane directions. We note there that a previous 
knowledge (MEM observation in that case) of the magnetic distribution could be a key 
point for such a quantitative analysis. 

The scattergram (inset of FiglUC) exhibits two main characteristics : a round shape due 
to the distribution of isotropic magnetic moments, and a center of shape which does not 
correspond to the center of the scattergram image. When re-centering the distribution 
of magnetic vectors in the scattergram image, the electrical contribution to the phase is 
removed (FigHlD). The small phase variations due to magnetic deflections are now easy 
to observe. 



7 



Another simple way of removing this component was previously proposed by Volkov 
and Zhu in |19|]. The aim of this method was to estimate the electrostatic phase by fitting 
a parameter corresponding to the log-ratio equation |23l] . This method also needed only 
one phase information and no previous quantitative thickness datas. Nevertheless, it 
is necessary to add another mathematical approach, namely fitting the parameter for 
electrostatical phase shift approxiamation. Furthermore this technique is only working 
with sharp edges in the imaging area. It could not have been implemented here. The 
two techniques discussed here can thus be described as complementary techniques for 
magnetic quantitative analysis in LTEM. 

We now discuss about the MIP variations in the image. The large field of view used 
in LTEM observation (around 1 [^m in that case) and the geometry of the TEM sample 
(prepared with traditional ion milling) leads to smooth variations of the MIP in the 
image. As a simple evaluation, the MIP here could not be evaluated with a precision 
of more than 1 V (see for example literature review for MIP of semiconductors in ref. 



251) . Over the field of view of that observation the MIP is mainly modified by the ratio 



of MgO over the global MIP of the magnetic layer. Considering the 6 degree angle 
used for sample preparation, the variation in the MIP value corresponds to less than 
1 V (considering MIPm^o = 14 F and MIFpePd = 25 V) which is less than the precision 
we could expect on such a MIP estimation. Then small variations of the MIP (due to 
grain boundaries, atomic steps or structural defects) are not considered in our experiment. 



We have then determined the pure magnetic induction map. A zoomed area of the 
reconstruction is shown in Fig|3l In that representation we observe two main areas where 
the magnetization is in-plane : 

• "Magnetic tubes" which are running all around the sample. These are the domain 
walls (intense colour on the map Fig|3l, enhanced by long thin arrows). 

• Other areas where the magnetization is also in-plane with a modulus which is 
varyng along the layer (and located around the domain walls) correspond to the 
domains (weaker colour on the map Fig|3l enhanced by small, wide arrows). The 
signal does not come from the out-of-plane domains but in the soft layer where the 
longitudinal magnetization area is parallel or antiparallel to the domain walls (mark 
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3 on Fig HI). 

We might expect some dark areas corresponding to pure perpendicular magnetization. 
These areas are not localised in the middle of the domains (due to the soft layer) but 
on the boundaries between domains and domain walls (see dashed areas in Fig [T]). The 
no-deflection condition can be explained by a balancing between the Bloch domain wall 
(mark 1 on Fig[T]) and the longitudinal magnetization area under the domain (mark 3). 
This balancing occurs only if the domain wall is running in the opposite direction to this 
longitudinal magnetization area. The dark areas on the color map (Fig [S]) do not always 
surround domain walls but only when a domain wall and the longitudinal magnetization 
area under the domain near this wall are antiparallel. 

If we assume that stray fields and Neel Caps are cancelling one another on both sides of 
the sample, we can evaluate magnetization inside the domain walls. Thus, considering a 
constant magnetic thickness of 37 nm at the domain wall, we found a value of magnetic 
induction around 1.3 T in the domain walls (long thin arrows in Fig |3l). This value 
is according to the FePd(Llo) saturated magnetization (measured to be 1.29 T in |24l]). 
On the contrary, considering a magnetic thickness of 16 nm between the domain walls 
(the magnetic induction in the LIq layer is purely perpendicular) we measured a value 
betweenn 0.01 and 1.0 Tesla in the soft layer (small and wide arrows in Fig|3l). These 
values are in the range expected for the disordered FePdi saturated magnetization. They 
could be over-estimated when measured close to the domain walls. 
Domain wall widths appear to vary, depending on shape (straight or curved) and the 
proximity of other domain walls. We can assume a size around 10 nm for the domain 
walls in this sample. Theorical values of 8 nm ITi] were calculated using an infinite 
problem without considering stray field energy or the shape of the thin foil. Neverthe- 
less, this kind of value has to be evaluated for different values of defocus by means of 
zero-defocus approximation. MFM measurements previously showed a pseudo-period 
for this sample of 160 nm, which is the same as our measurements with an auto-correlation. 



Along the walls, "switching" points were observed where the magnetization abruptly 
turns 180° with respect to the wall direction. One of these points is indicated in the Fig 
|3] with a dashed circle. It can also be observed in the Fresnel contrast as an inversion 
from dark-bright contrast to bright-dark contrast. These singularities could be of great 

9 



interest if it could be confirmed that they are Vertical Bloch Lines |26l, 1270 . The Bloch 
lines observed in traditional in-plane anisotropy materials play a major role during the 
magnetization process. This role will be discussed in a future publication on the nature 
and the behaviour of these lines. 



Conclusion 



We have demonstrated that Lorentz Microscopy is potentially a powerful tool for the 
investigattion of magnetization configuration at the magnetic domain walls scale. It is 
possible to study materials which are anisotropic out-of -plane without tilting the speci- 
men. FePd alloys have been studied and their Bloch wall structure revealed. We have 
noted some defects that appear along the domain walls which are potentially impor- 
tant during the magnetization process. Observations have been carried out using a new 
method for subtracting mean inner potential information from phase datas. This subtrac- 
tion enables the quantification of the magnetic map with a good correlation to previous 
measurements. Resolution of the Transport of Intensity Equation is a very high-resolution 
technique both spatially (around ten nanometers) and in terms of magnetic sensitivity 
(few tens of nm.T). 
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Appendix 



We describe here the notion of scattergram which can be labeled as a bidimensionnal 
histogram of the moment vectors by gray level scale. As shown in Figure |5l we use a 2D 
representation to count the magnetic moments in a vectorial map. Starting with two 
images, displaying two components of magnetization (see Fig |5] upper left and upper 
right) a vectorial representation can be made (see Fig |5] bottom left). Each pixel in the 
scattergram (See Fig|5]bottom right) represents a couple of X-intensity and Y-intensity but 
also a couple of modulus and angle. The origin (null modulus) is choosen as the center 
of the image. Finally the intensity of each pixel of the scattergram represents the number 
of times the corresponding couple {{modulus, angle} or {intensityX, intensityY}) appears in 
the map. If we sum all the values of the scattergram image, we get the number of pixels 
used for the inital images (mostly 1024x1024 = 1048576). 
Filling of the scattergram image can be described by the following equation : 

S(z, ;■) = P {z < m^ < z + 1 & < my < + 1} (8) 

where the S(z, ;) are the pixels of the scattergram and the (m^, my) are the components of 
magnetization normalized by the pixel size. P is the probability of a {m^, my) occurence 
in the image. The scattergram can be also seen as a representation of the back focal plane 
of the Lorentz lens, the origin of the scattergram being the transmitted beam. Interested 



readers can also find information in |28|1 
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Figures 




Figure 1: Micromagnetic simulation from GL_FFT on the FePd bi-layer. The input parameters are 
displayed in the text. The simulations contain only the magnetic layers. The growth direction is 
along the vertical. A colour scale is used for the component of magnetization perpendicular to the 
figure. Marks on the figure show the in-plane magnetic components described in the text. Dashed 
areas correspond to some areas which fulfil the no-deflection conditions. 
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Figure 2: A. In-focus image of FePd sample. B. Over-focus (ISfim) of the same area as in A. The 
inset is a magnified view of the contrast of a domain wall. Arrows pointing the domains were 
choosen arbitrarily. C. Idealized profile of the perpendicular magnetic anisotropy system in plane 
view geometry. Domains are separated by Bloch walls. The marks 1-3 correspond to the same 
magnetic components as in Fig[TJ The same color scale as in Fig[l]is also used. B. Corresponding 
Fresnel contrast obtained by simulation. 
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Figure 3: Magnetization map from a phase retrieval process issued from a zoom in Fig HI The 
color scale used here shows the magnetic vector direction with a colour and the intensity of the 
colour corresponds to the modulus of this magnetic vector (see scale as inset). Thin and long 
white arrows are pointing the magnetization in the domain walls. Thick and small white arrows 
are pointing the magnetization laying under the domains. The dashed circle points a singularity 
in the wall described in the text. 
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Figure 4; A. Under-focused (-20 [im) image used for the phase retrieval process. B. Corresponding 
over-focused (20 [im) image. C. Phase image obtained after the reconstruction from focal serie. 
The inset is the scattergram associated to the phase gradient (in temperature scale). D. Phase image 
and its corresponding scattergram after electrical contribution removal. E. Explicative scheme of 
the phase shift origin. The phase of the electron wave (red line) is shifted by an electrical and a 
magnetic origin. The decomposition of (p and its gradient (directly related to B) is explained. F, 
Fresnel contrast simulation (A/ = 20 /^im) obtained from the magnetic induction map determined 
from phase image. The domain walls appear as coupled dark-bright bands. 
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Figure 5: Scattergram for a 16-pixels image representing a schematic Cross-Tie wall. The two top 
images are the X and Y components of the magnetization. The bottom left image is a vectorial 
image representation of the magnetization distribution.The bottom right image is the scattergram. 
The color scale used in the scattergram is zero for black and the gray levels correspond to the 
among of times the corresponding vector is dispalyed in the previous images. Colors are used to 
show the position of the vectors in the scattergram. 
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